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Abstract 

An antiferromagnetic version of the well-known majority voter model on square 
and honeycomb lattices is proposed. Monte Carlo simulations give evidence 
for a continuous order-disorder phase transition in the stationary state in both 
cases. Precise estimates of the critical point are found from the combination 
of three cumulants, and our results are in good agreement with the reported 
values of the equivalent ferromagnetic systems. The critical exponents l/v, ^ jv 
and Pjv were found. Their values indicate that the stationary state of the 
antiferromagnetic majority voter model belongs to the Ising model universality 
class. 
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1. Introduction 

The Majority voter model (MVM) is a simple non-equilibrium Ising-like 
system, proposed as a way to simulate opinion dynamics. The collective be¬ 
haviour of the voters shares many aspects with the well-established theory of 
non-equilibrium phase transitions and results form simulations can be analysed 
similarly [l[. In the standard MVM, the system evolves following a dynam¬ 
ics where each “voter” (spin) assumes the same opinion as the majority of its 
neighbours, with probability (1 -I- x)/2 and assumes the opposite opinion, with 
probability (1 — x)/2. Here x is the control parameter, with a range 0 < x < 1. 
The stationary state of the MVM presents a second-order phase transition, at 
some critical value Xc, and previous numerical works on regular lattices show 
that the critical exponents are compatible with its equivalent lattices for the 
Ising model 0, S i 0. Those results seem to confirm the conjecture that 
non-equilibrium models with up-down symmetry and spin-flip dynamics fall in 
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the universality class of the Ising model However, numerical simulations in 
non-regular latices, such Archimedean, small-world or random lattices, rather 
seem to indicate that the critical exponents are governed by the lattice topol¬ 
ogy EiiE 


111 112 , 1131 , 114 , 1151 . 


On the other hand, the MVM belongs to a family of generalized spin mod¬ 
els la] that can be modelled in terms of a competing dynamics, induced by heat 
baths at two different temperatures (on two dimensional square lattices) [iZl, llSj 
and hence have a non-equilibrium stationary state. The equilibrium stationary 
state of the ferromagnetic Ising model is recovered when both temperatures 
are equal; in this case, the detailed-balance condition is fulfilled. All members 
of the family share the same definition of the order-parameter and there is a 
critical line that separates the paramagnetic (disordered) phase from the ferro¬ 
magnetic (ordered) phase 0. So far, the ferromagnetic-paramagnetic transition 
out of equilibrium has been extensively studied with a single-flip dynamic rule 
that recovers the equilibrium Ising model. When antiferromagnetic interactions 
are included in the Ising model, the phenomenology becomes richer, additional 
phases with different critical behaviour, multicritical points, etc. can appear (see 
Ref. 0 and references therein). In this work, we want to explore, as a start¬ 
ing point, if it is possible to implement an antiferromagnetic non-equilibrium 
version that follows a single-flip updating scheme. Godoy and Figureido 0 
proposed a non-equilibrium mixed-spin antiferromagnetic model with two up¬ 
dating schemes: one single-spin via Glauber dynamics and a two-spin updating 
linked to an external energy source. The aim of the present work is to implement 
the antiferromagnetic version of the MVM on the square and the honeycomb 
two-dimensional lattices, evaluate the critical point and the stationary critical 
exponents, 7 and v. 

This work is organised as follows: in section 2, the antiferromagnetic MVM 
is defined and the hnite-size scaling method used to analyse its stationary state 
is outlined. In section 3, the results of the Monte Carlo simulation are reported 
and the critical parameters are extracted. We conclude in section 4. 


2. Model and Finite-Size Scaling 

In the ferromagnetic version of the MVM Q , each lattice site is occupied by 
a spin, (Ti, that interacts with its nearest neighbours. The system evolves in the 
following way: during an elementary time step, a spin ( 7 ^ = ±1 on the lattice is 
randomly selected, and flipped with a probability given by 


p{x) 


-{l + px). 


( 1 ) 


50 Herein, 77 stands for the (ferromagnetic) rule 


f -sgn[i7j • cTi] ; it ^ 0 
\ 0 ; if iJ. = 0 


and where Hi is the local field produced by the nearest neighbours to the i th 
spin. The control parameter x acts analogously to a noise in the system. With 
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the dynamics (fim . a given spin di adopts the sign of Hi (the majority of 
its nearest neighbours) with probability (1 + x)/2 and the opposite sign of Hi 
(the minority) with probability {1 — x)/2. In the Ising ferromagnet, the sign 
of the bilinear exchange interaction in the hamiltonian defines if the system is 
ferromagnetic, or antiferromagnetic, respectively, for a positive or negative sign. 

An ontiferromagnetic version of the MVM is now naturally obtained by 
replacing the rule ([2|) by 


J +sgn[H, ■ di] ; it H,^0 , , 

^ \ 0 ; if iJ, = 0 ^ ’ 

in combination with the flip probability o, with 0 < a; < 1. The anti-ferro 
MVM is defined by the dynamics (fira . with the control parameter 0 < a; < 1. 
Alternatively, one could keep the rule but replace x i—>• —x in the flip 
probability (HD. The consideration of square and honeycomb lattices allows to 
study the effect of having an even or odd number of nearest neighbours of each 
spin on the critical behaviour in the stationary state. 

In analogy with simple Ising magnets, in order to measure the paramagnetic- 
antiferromagnetic phase transition, we shall use the staggered magnetisation as 
order-parameter 



where N is the total number of lattice sites, Ci takes values of ±1 depending in 
which sub-lattice the site is located (see Fig.[T|) and (...) stands for time average 
taken in the stationary regime. 


a) 


-t—V 


-V —f- 


b) 



Figure 1: Sub-lattices for a) the square and b) the honeycomb geometries. Here, Cj = 1 
for white sites and = — 1 for black sites. In the square lattice we use periodic boundary 
conditions and in the honeycomb lattice we use use skew boundary conditions in the horizontal 
boundaries. 
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The (staggered) susceptibility is given by 


X = Nx{{m?) — {m)^}. 


( 5 ) 


We shall use the method proposed in Ref. [2l|, where three different cumulants 
are used for the evaluation of the critical point: (i) the fourth-order or Binder 
cumulant 221 

( 6 ) 


= 1 - 


3(m2)2^ 


(ii) the third-order cumulant (where (m^) is defined analogously to eq. (|4])) 


{ 7 ( 3 ) = 1 


{w?) 

2 (m)(m 2 ) ’ 


and (iii) the second-order cumulant 

It{2) _ ^ 2(m^) ^ 


( 7 ) 

( 8 ) 


The scaling forms for the thermodynamic observables, in the stationary state, 
and together with the leading finite-size correction exponent uj are given by 


m(e, L) s 

s L-^/''{M{eL^/'') + L-^M{eL)), 

( 9 ) 

X(C A) s 

s L^/‘^(x(eA'/‘^) + A-"^(eA)), 

( 10 ) 

C/(P)(cA) s 

^ [p] 

s U^P\eL^/'')+L-^U (eL). 

( 11 ) 


where e = x—Xc is the distance from criticality, p = 2, 3 or 4, and L = '/N is the 
linear size of the lattice. The parameters /3, 7 and v are the critical exponents 
for the infinite system, see [I| for details. 

In principle, the critical point Xc is found from the crossing points in the 
cumulants A precise estimation of Xc is achieved by taking into account 

the crossing points for different cumulants and with p ^ q arise for 
different values of L. The values of x, where the cumulant curves U^^\x) for 
two different linear sizes Li and Lj intercept are denoted as x[^'^. We expand 
Eq. (ITT]) around e = 0 to obtain 


uip) ~ t/b) + U^p)(:L^/'' + -7 0(e2, eL-^) 


( 12 ) 


where are universal quantities, but [7^^^ and are non-universal. The 
value of e where the cumulant curves for two different linear sizes Li and Lj 
intercept is denoted as e^fj . At this crossing point the following relation must 
be satisfied: 

(13) 




Here := U /U^p\ Combining for different cumulants {q ^ p) we get 


Ap) I Aq) 

b' - T - 4 

— Xc \X^, 


(14) 


4 







where Apq = and is non-universal (see Refs. 21 


1^ for additional details). Equation (IT^ is a linear equation that makes no 
reference to or w and requires as inputs only the numerically measurable 
crossing couplings y The intercept with the ordinate gives the critical point 
location. 


3. Results 

100 As a first step, we check that the dynamics Eqs. (ITI31) gives rise to a pa- 
ramagnetic-antiferromagnetic phase transitions, depending on the value of the 
parameter x > 0. For a qualitative illustration, in Fig. [5] we present snapshots 
for two different values of x, of the state of a square lattice of size L = 256. 
Clearly, the two states look very different and are also distinguished by the 
respective values of (m) as defined in (g]); this suggests the existence of a phase 
transition, at some intermediate critical point Xc- 



Figure 2: Snapshots of the square lattice with size L = 256, for a. x = 0.1 (left) and x = 0.95 
(right). The order-parameter values are (m) ~ 0.0034 and (m) ~ 0.9446, respectively. 


We performed simulations on three different lattices with linear sizes L = 24, 
28, 32, 36, 40 and 48; and did this for both the square and honeycomb lattices, 
see figure [T] Starting with a random configuration of spins, the system evolves 
following the dynamics given by eqs. ma . Then, we let the system evolve during 
a transient time, that varied from 4 x 10^ Monte Carlo time steps (MCTS) for 
L = 24 to 1.2 X 10® MCTS for L = 48. Averages of the observables were 
taken over 2 x 10® MCTS for L = 24 and up to 1.2 x 10^ MCTS for L = 48. 
Additionally, for each value of x and L, we performed up to 400 independent 
runs, in order to improve the statistics. 

In Fig. [3] we show the third-order cumulant curves as function of x for the 
square and for the honeycomb lattice. The hgure shows clearly that in both ge- 
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ometries the curves for different sizes cross around a certain (lattice-dependent) 
value of X. Similar behaviour has been observed as well for and for 



X 


Figure 3: Third-order cumulant as function of the parameter x in (a) the square and (b) 
honeycomb lattices. The solid lines are the third-order polynomial fits. 


For the evaluation of the critical points, we used a third-order polynomial 
fit for the cumulant curves. Recalling eq. da, the estimation of the critical 
points for the two cases are shown in Figure 01 where we plot the variable 
O' := over against the variable 5 := We can observe 

that the raw differences between crosses, S, are smaller in the square case. But, 
since the smaller values in S are around for the square lattice and for 
the honeycomb lattice (these values correspond to the crossings between the 
largest sizes using in our simulation), we can be sure that largest sizes will not 
improve significantly our results. In any case, the largest source of numerical 
error apparently comes from the statistical uncertainty of the data points. If 
we take another pair of cumulants, the dehnitions of a and 6 are adapted in 
an obvious way. The linear fits of Eq. (ITT)) give the following estimates for the 
critical points 


Xc = 0.84969(4) 
Xc = 0.87195(14) 


6 


; square lattice 
; honeycomb lattice 


(15) 
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0.8497 


0.8496 
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Figure 4: (Color online) Evaluation of the critical point Xc for (a) the square lattice and (b) 
the honeycomb lattice. The points represent the numerical data obtained from third-order 
polynomial fits and the dashed lines are the linear fits of Eq. m. Smaller <5-values correspond 
to larger system sizes. Black circles: crossing of red squares: crossing of 

blue diamonds: crossing of 


where the numbers in brackets give the estimated uncertainty in the last given 
digit (s). Both results are in good agreement with the reported values for the 
ferromagnetic MVM We want to point out that, when we compare the 

data for the crossing of the curves with the previous reported data 

for the ferromagnetic case from Ref. the range in the differences <5 is almost 
two times larger in the ferromagnetic case. Since the linear sizes considered are 
the same in both cases, we suspect that the scaling effects are smaller in the 
antiferromagnetic case. Additional simulations in the Ising model, for antiferro¬ 
magnetic and ferromagnetic, interactions for different lattice geometries and for 
the MVM on square lattices would help to see if this is a universal behaviour. 

The critical exponents can be evaluated, by using Eqs. (HU, at the critical 
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point e = 0. One expects the following finite-size scaling behaviour 


m{L) 

oc 

(16) 

X{L) 

oc 

(17) 

dU^P'i 

oc 

X — Xc 

(18) 

dx 


In Fig. [5l we show the derivatives of the cumulants at the critical point. From 
the Hnite-size scaling law (flSl) . we obtain the following results: Ijv = 1.02(4) for 
the three cumulants in the square lattice and l/iz = 1.02(4), 1.03(4) and 1.03(4) 
for U^‘^\ and respectively in the honeycomb case. The evaluation of 
150 jv is shown in Fig. [ 6 l with the relation (flTll we obtain jv = 1.756(8) and 
^jv = 1.759(10) for the square and honeycomb lattices respectively. We present 
in Fig. [7] the evaluation of /3/jz, with Eq. (fT6l) we obtain fijv = 0.123(2) and 
^jv = 0.124(8) for the square and honeycomb lattices, respectively. 

Our results are collected in table [H where we also included the correspond¬ 
ing results for the ferromagnetic MVM ii- Clearly, all reported numerical 
estimates are very close to the exactly known values of the two-dimensional 
Ising ferromagnet, e.g. [l|, appendix A]. According to hyperscaling, one expects 
(2,5 -b 7 )/j^ = d = 2, which is well conhrmed by the simulation results. We also 
see that the agreement for the anti-ferromagnetic MVM is even slightly better 
than for the ferromagnetic MVM. It turned out that an explicit consideration 
of possible finite-size corrections, described by the Wegner exponent w, was not 
necessary. 


Table 1: Values of the critical parameters for the antiferromagnetic Majority voter model 
(AMVM). For comparison, we also include the values reported for the ferromagnetic Majority 
voter model (FMVM) by Kwak ei al. [3| on square lattices and by Acuha-Lara et ai Q on 
honeycomb lattices. 


Xc 

l/v 

llv 

filv 

(7 + 2/3)/z/ 

model 

0.84969(4) 

1.02(4) 

1.756(8) 

0.123(2) 

2.002(9) 

Square AMVM 

0.87195(14) 

1.03(4) 

1.759(10) 

0.124(8) 

2.007(15) 

Honeycomb AMVM 

0.8500(4) 

0.98(3) 

1.78(5) 

0.120(5) 

2.02(5) 

Square FMVM 

0.8721(1) 

1 .01(2) 

1.755(8) 

0.123(2) 

2.001(9) 

Honeycomb FMVM 


Additionally, we evaluated the universal quantities (lattice-dependent) 
for both lattices, using Eq. m with e = 0. Our data do not allow to reliably 
evaluate the Wegner exponent w, since there is no need to include scaling cor¬ 
rections in these systems, as we illustrate in Fig. El Our result for the Binder 
cumulant, = 0.612(1), for the square lattice is in good agreement with the 
reported values of the two-dimensional Ising model, = 0.610690(2) [ 2 ^, or 
= 0.6106922(16) [ 2 ^. We listed all the results for the three cumulants, for 
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Figure 5: (Color online) Log-log plot of the derivatives of the cumulant (a) Idx, (b) 

gf/(3) /partialx and (c) /partialx, taken at the critical point x = Xc- Black circles 

correspond to the square lattice and red squares to the honeycomb lattice. The dashed lines 
are power-law fits. 


both lattices, in Table [H We observe that the results for the honeycomb and 
square lattices are clearly different from each other, in all cases. The expla¬ 
nation for this behaviour is simple: because in our simulations, the boundary 
conditions of the two lattices are different: we use periodic boundary conditions 
for the square lattice and skew conditions in the horizontal boundaries for the 
honeycomb lattice (see Fig.[T|). Estimates of reported in the literature, and 
based on simulations in two-dimensional critical Ising model, produced a simi¬ 
lar spread for different lattices or boundary conditions [1^ . In order to test the 
universality through the cumulants it would be necessary to perform sim¬ 
ulations for the Ising model on the honeycomb lattice with the same boundary 
conditions. 
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Figure 6: (Color online) Log-log plot of the susceptibility at the critical point for the square 
(black circles) and honeycomb (red squares) lattices. The dashed line are power-law fits. 


4. Conclusions 

We have introduced an antiferromagnetic version of the MVM, both on the 
two-dimensional square and honeycomb lattices, and studied its stationary state, 
as a function of the control parameter x, through intensive Monte Carlo simu¬ 
lations. The phase transition in the stationary state, between an ordered phase 
for X large enough and a disordered phase for x small enough, falls on both lat¬ 
tices into the 2D Ising model universality class, and independently of whether 
the number of nearest neighbours is even or odd. This model further illus¬ 
trates that the set of critical exponents is consistent in two-dimensional regular 
lattices. Future work should explore if additional phases and/or multicritical 
points exist if an external field is included. 
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Figure 7: (Color online) Log-log plot of the order parameter at the critical point for the square 
(black circles) and honeycomb (red squares) lattices. The dashed lines are power-law fits. 
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